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Abstract. Motivated by experimental results on we use a factorization approach to study 

these decays. Two mechanisms concerning kaon pair production arise: current-produced (from vacuum) 
and transition (from the B meson). The kaon pair in the B'^ — > D^*^^ K~ K'^ decays can be produced 
only by the vector current (current-produced), whose matrix element can be extracted from e'^e~ KK 
processes via isospin relations. The decay rates obtained this way are in good agreement with experiment. 
The B~ D^*^'^ K~ K'^ decays involve both current-produced and transition processes. By using QCD 
counting rules and the measured B~ D^*^^^ K~ K^^ decay rates, the measured decay spectra can be 
understood. 
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1 Introduction 

The B K^*^^ decays have been observed for the 

first time by the Belle Collaboration T, with branching 
fractions at the level of lO""* — lO^'^. Angular analysis re- 
veals that K~K^ and K~ K*^ are dominantly = 1~ 
and 1+, respectively. While there is no sign of decay via 
resonance for the K^K^ pair, data suggest a dominant 
ai(1260) resonance contribution in the production of the 
K*'^ . The K~ K^*^'^ mass spectra are peaking near 
threshold. 

The near-threshold peaking of the K~K^*^^ mass spec- 
tra suggests a quasi two-body process where the colinear 
K~K^*^^ recoil against the D*^*-* meson. This is suggestive 
to apply fatorization to the three-body case j^j . Two kinds 
of decay amplitudes arise due to the flavor structures of 
the mesons: 0^*^+ R- K'-*^° involves {K-K'-*^°\V - 
A\0), with K~K^*'>^ produced by a weak V — A current; 
£)(*)o^-;^(*)o involves {K-K^*'^°\V - A\B-), where 
goes into K~K^*'>'^ via a weak current. 

In {K-K°\V-A\0), the K'K^ can only be produced 
by the vector current, and should be dominantly 1~. By 
isospin rotation, the kaon weak form factor {K~ K^\V\0) 
can therefore be related to the kaon electromagnetic (EM) 
form factors in e+e^ annihilation, where much data ex- 
ist. One can then calculate the rate without any tuning 
parameters. The predicted K~K° mass spectrum can be 
shown to have a peak near threshold, which arises from 
the kaon form factor and can be checked by experiment. 

In B- D(*^"K-K° decays, the R-R^' can also be 
produced by a current that induces B~ R^ R'^ transi- 



tions. The relevant matrix element {K~R^\V — A\B~) 
is parameterized by several unknown form factors, due 
to which the rate cannot be calculated. Nevertheless, by 
using a naive parametrization based on QCD counting 
rules (3j, the parameters in these unknown form factors 
can be determined from the decay rates, and the decay 
spectra can be obtained, which again have threshold en- 
hancement as closely related to QCD counting rules, and 
can be tested experimentally. 

The R-R*° in W ^ D'^*^+R-R*° can only be pro- 
duced by a weak current. The experimental observation 
that R~ R*'^ is in l"*" suggests a dominant axial current 
contribution. Although the decay rate cannot be calcu- 
lated due to the absence of data of the R~R*'^ axial form 
factors, it has been proposed that one can extract the 
j{-j^*o g^xial form factors given the R~R*° mass spec- 
trum 

In what follows, wc shall concentrate on decays that 
involve only R^ R'^ . In the next section we will introduce 
the relevant formalism and describe the numerical results. 
A discussion will be given in the last section where the 
conclusion is drawn. 



2 Factorization Formalism 

Starting with the relevant effective Hamiltonian, and the 
factorization ansatz, one arrives at |2| 

A{D^*^+R-R^*^^) = ^V,,V:,ar{D'^*^+\{cb)v-Am 
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Table 1. B{B° D'^*^+K-K°) in units of 10"''. 



MS BSW Experiment Q] 

W^D+K-K° l-&7lolt l-54lo2o 1.6 ±0.8 ±0.3 
_ ' ' < 3.1 (90% CL) 

B°^D*+K~K° 2.8t°Qll 3-05lo 39 2.0 ±1.5 ±0.4 

< 4.7 (90% CL) 
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Fig. 1. Fit to timelike 1-?^+ I (upper) and l-F^ol (lower) form 
factor data, where the inset is for the (j> region. 

x(X-xMO|(Ju)y|0), (1) 

A{D(*'>"K-K^*'>") = ^V,hV:,[a,{D(*'>"\{cb)v-A\B-) 
v2 

x(K-K''*^°\{du)v\0) 
+a2{K-K^*^\db)v-A\B-) 
x{D'^*^'>\{cu)v-Am (2) 

in which 

{K-{pK-)K^{pKo) \V^\ 0) = [PK- -PKoYF^^iq^) (3) 

in the isospin hmit, and 

{K- {pK-)K\pKo)\{V - A)^\B-{pb)) 
^ iw-{q^){pK- -PK«)t, 

+ Hq^)e^^a.pp''Bq°'{pK- -pKof, (4) 

where q = px- + Pa'"- The fact that is in 1^ has 

been taken into account in the above parametrizations. 
The {D^*'^\{ch)v-A\B) is the same as in two-body cases 
and we adopt both the BSW and the MS [H] models 
for comparison. 

The kaon weak vector form factor is related to 

its EM partners via the isospin relation 

F^''{q^) = FK+{q^)~FKo{q^), (5) 

where , F^o are the EM form factors of the charged 
and neutral kaons, respectively. By fitting to the EM data, 
one can obtain the kaon EM form factors and hence the 



weak vector form factor |2], as shown in Figs. ^ and 13 
Readers are referred to [2] for detailed discussion on the 
fitting of kaon EM form factors. 

For the unknown form factors in Eq. 0] we take the 
following naive parametrization [5j 

= M0 = |, (6) 

where c^ji are free parameters to be fitted by data. The 
1 /t^ arises from the minimum number of hard gluons to 
produce an energetic kaon pair from a decaying B meson, 
which characterizes the asymptotic behavior of the form 
factor. 

In Tabled we show the calculated branching fractions 
of the W D^*'^+K-K^ modes. One can see that the 
results are in very good agreement with experiment. In 
the absence of any tuning parameters in our formalism 
for these decay modes, the agreement between experiment 
and the model actually provides evidence for factorization! 

On the other hand, the predicted K~ mass spec- 
tra in Fig. 13 of the D^*^^ K~ K'^ modes show peaks close 
to threshold, which is due to the near-threshold behavior 
of the F{^^ form factor (see Fig. EJ. The re is no other 
clear structure, other than the B D* form factor ef- 
fect at larger q^ . Because of lower £)(*)+ reconstruction 
efficiencies, the spectra has yet to be measured experi- 
mentally PP, but our predicted spectrum can be checked 
soon with more data. 

For B~ D^*^'^K~ K'^ which involve the unknown 
parameters Cw_, Ch, we fit to the decay rates to obtain 
their values, as shown in Table El The K~K^ mass spec- 
tra can then be obtained, in particular for the B~ — > 
D^K~K^ mode, where one can see from Fig. 0] that the 
K~ mass spectrum can be roughly described by the 
naive form factor model of Eq. The mass spectrum of 
B~ D*^K~ in Fig.Elcan be checked in experiment. 
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Table 2. Fitted values of transition form factor parameters c-u,_ and Ch, by using the central values of D'^*^^'' K" K'^ rates. See 
[5] for detailed discussions. 



MS(BSW) 

C-w_ 


(GeV') 


MS(BSW) 


(GeV') 




B- D°K-K° -35.4 (- 
B- D*°K-K° -35.4 (- 


-33.0) 
-33.6) 


11.3 (13.1) or - 


-16.1 (-18.5) 


5.5 ±1.4 ±0.8 
5.2 ± 2.7 ± 1.2 




Fig. 3. The K' K° mass spectrum for B° D+ K' K° 
(lower) and D*'^ K~ (upper), where solid (dashed) line 
stands for using the MS (BSW) hadronic form factors. 
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Fig. 4. B D^K spectrum, where solid (dashed) line 
is for the MS (BSW) model, and the data is from Ref. [T]. 
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Fig. 5. The K' K° mass spectrum for B' K' K° , 

where solid, dot-dashed, dashed and dotted lines are for MS 
model with cu = 11.3, -16.1 GeV^ and BSW model with 
Ch = 13.1, —18.5 GeV^, respectively. 



3 Discussion and Conclusion 

A factorization approach has been used to study three- 
body S JsT" JsT" decays. There are two mechanisms 
of kaon pair production, namely current-produced and 

transition. For B — > D^*^~^K~ which involve the for- 
mer, one can make use of kaon EM data through isospin 
rotations. The result is in good agreement with experi- 
ment, which supports factorization. 

The — > D'-*^^K~ decays also receive the transi- 
tion contribution. The form of these transition form fac- 
tors are determined through QCD counting rules, and we 
fix the parameters by using the measured D^K~K^ and 
D*^K~K^ decay rates. The predicted mass spectra of the 
D^K^K^ mode agrees well with data and ex- 
hibit threshold enhancement as do the B° D^*^^K~K^ 
cases. Despite the success in describing the mass spectrum 
of the D°K-K° mode, our treatment of the B- R-R^ 
transition form factors may be oversimplified. Assuming 
the asymptotic form required by PQCD may be too strong 
an assumption, and might have over-enhanced the contri- 
bution from the near-threshold region. More careful study 
on other possibilities, such as using pole models for tran- 
sition via resonances, would be helpful in clarifying the 
underlying dynamics of the B^ — ^ D'^*^^R^ R^ transi- 
tions. 

Finally, with the feasibility of extracting R~ R*'^ axial 
form factors emboldened by the success in B -> D'^*^R-R° 
decays, B decay data plus factorization have opened up a 
new avenue to the study of meson form factors, which have 
traditionally been fundamental quantities to many fields 
in both nuclear and elementary particle physics. The suc- 
cess of factorization in B D^*^ R^ R^ decays urges a 
serious study of the underlying mechanism. 



References 

1. A. Drutskoy et al. [Belle Collaboration], Phys. Lett. B 542, 
(2002) 171 arXiv:hep-ex/02 07041, . 

2. C. K. Chua, W. S. Hon, S. Y^iau and S. Y. T sai, Phys. 
Rev. D 67, (2003) 034012 arXiv:hep-ph/0209164l. 

3. S. J. Brodsky and G. R. Farrar, Phys. Rev. D 11, (1975) 
1309. 

4. C. K. Chua, W. S. Hon, S. Y. Shiau and S. Y. Tsai, in 
preparation. 

5. M. Wirbel, B. Stech and M. Bauer, Z. Phys. C 29, (1985) 
637. 

6. D. Melikhov and B. Stech, Phys. Rev. D 62, (2000) 014006 
[arXiv:hep-ph/0001113 . 



